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AN ABSTRACT ANALYSIS OF 
OPTIMAL GOAL-ORIENTED ADAPTIVITY 

MICHAEL FEISCHL, DIRK PRAETORIUS, AND KRISTOEEER GEORGE VAN DER ZEE 


Abstract. We provide an abstract framework for optimal goal-oriented adaptivity for 
finite element methods and boundary element methods in the spirit of [15]. We prove 
that this framework covers standard discretizations of general second-order linear elliptic 
PDEs and hence generalizes available results [7, 37] beyond the Poisson equation. 


1. Introduction 

1.1. State of the art contributions. This work considers the simultaneous adap¬ 
tive control of two error estimators and which satisfy certain abstract axioms 
from Section 2.4, below. The estimator product is designed to control the error in 

goal-oriented adaptivity and allows to prove optimal error decay for the goal functional. 
This is discussed in Section 1.2 and demonstrated in Section 4-6 for various model prob¬ 
lems. We analyze three adaptive mesh-refinement algorithms (Algorithm A, B, C) which 
allow optimal convergence rates for the estimator product in the sense that each possible 
algebraic convergence r > 0 will be achieved, i.e., 

Vu,iVz,i < (#7^ - for all i e N, 

without any a priori knowledge of the optimal rate Topt- Here, the triangulations {Ti)e£No 
are generated by the respective adaptive algorithm starting from some given initial tri¬ 
angulation To- 

While standard adaptivity aims to approximate some unknown exact solution u at op¬ 
timal rate in the energy norm (see, e.g., [18, 23, 42] for adaptive FEM, [21, 22, 24, 26] for 
adaptive BEM, and [15] for a recent overview on available results), goal-oriented adap¬ 
tivity aims to approximate, at optimal rate, only the functional value g{u) (also called 
quantity of interest in the literature). Goal-oriented adaptivity is usually more important 
in practice than standard adaptivity. It has therefore attracted much interest also in the 
mathematical literature; see, e.g., [6, 8, 9, 19, 28, 31, 40] for some prominent contribu¬ 
tions. However, as far as convergence and quasi-optimality of goal-oriented adaptivity 
is concerned, earlier results are only the two works [7, 37] which are concerned with the 
Poisson model problem and the work [29] which considers general second-order linear el¬ 
liptic PDEs, but is concerned with convergence only. Moreover, the analytical arguments 
of [7, 37] are tailored to the Poisson equation and do not directly transfer to the more 
general setting of [29]. The quasi-optimality analysis for goal-oriented adaptivity has also 
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been named as an open problem in the recent work [13]. In view of this, the contributions 
and advances of the present work can be summarized as follows: 

• We give an abstract analysis for optimal goal-oriented adaptivity which applies 
to general (non-symmetric) second-order linear elliptic PDEs in the spirit of [23] 
which even extends the problem class of [29]. 

• The analysis avoids any (discrete) efficiency estimate and thus allows for simple 
newest vertex bisection, while [7, 37] require local bisec5-refinement in the spirit 
of [42]. 

• Unlike [7], our proofs avoid any assumption on the resolution of the given data 
as, e.g., a saturation assumption [7, eq. (4.4)]. 

• Unlike [29], our analysis does not enforce the condition that the initial triangu¬ 
lation To is sufficiently small, since we do not exploit the regularity of the dual 
solution to prove some crucial quasi-Galerkin orthogonality. 

• Finally, our analysis does not only cover the finite element method (FEM), but 
also applies to the boundary element method (BEM). 

Related recent work includes [20], where the goal are point errors in symmetric BEM com¬ 
putations. Although we shall verify the mentioned estimator axioms only for standard 
FEM and BEM discretizations, we expect that they can also be verified for discretiza¬ 
tions in the frame of isogeometric analysis; see, e.g., [34] for some goal-oriented adaptive 
IGAFEM. 

1.2. Goal-oriented adaptivity in the frame of the Lax-Milgram lemma. The 

following introduction covers the main application of the abstract theory, we have in 
mind. Let A be a Hilbert space with norm || • \\x, and let o(-, •) : A x A —>■ M be a 
continuous and elliptic bilinear form on A. For a given linear and continuous functional 
/ G A*, let u G A be the unique solution of 

(1) a{u , v) — f{v) for all n G A. 

Let G A* be the so-called goal functional, i.e., we aim to approximate g{u) at optimal 
rate. To this end, suppose that associated with each triangulation % of some problem 
related domain U C M'^, there is a finite dimensional subspace A* C A. Let U* G A* be 
the unique Galerkin approximation of u which solves 

(2) a(U* , U*) = /(U*) for all U* G A*. 

Furthermore, let .2 G A be the unique solution to the so-called dual problem 

(3) a(v , z) = g(v) for all n G A. 

Then, for any Z* G A*, it follows 

(4) \g{u) - g{U^)\ = \a{u - U^,z)\ = \a{u -U^,z- Z*)| < ||u - U^\\x ||^ - ^*|U- 

Here and throughout, < abbreviates < up to some generic multiplicative factor G > 0 
which is clear from the context, e.g., the hidden constant in (4) is the continuity bound 
of o(-, •). Suppose that we compute the unique Galerkin approximation Z* G A* of the 
dual solution z G Al, i.e., 


a(K,Z,)=g(K) for all U*G A*, 


( 5 ) 
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and that the Galerkin errors on the right-hand side of (4) can be controlled by computable 
a posteriori error estimators 


(6a) 

11“ - Ufllx < Vu,* ■= 

(y: ,„,.(T)f 

Te% 

(6b) 

Ijz - Z^Hx < rjz,* ■= 1 



Te% 


Under these assumptions, we are altogether led to 

(7) \9{u) - g{U^)\ <Vu,*Vz,*- 

Overall, we thus aim for some adaptive algorithm which drives the computable upper 
bound on the right-hand side of (7) to zero with optimal rate. 

Remark 1. Using the residual Ru,-k{v) := f{v) — a({7*,'y) G X*, one ean improve the 
hidden eonstants in (4) by 

\g{u) - g{U^)\ = \a{u -U^,z- Z*)| < sup \\z - Z^\\x =: \\Ru,*\\x* lk - Z^\\x. 

\\v\\x 

Some a posteriori error estimators (such as, e.g., estimators based on flux equilibra¬ 
tion [44]^ allow for the reliability estimate of the residual ||Ru,*||a’» < C'rei?7u,* even with 
known constant Oei = 1- The same arguments for the residual RzAv) := q(v)—a(v, ZX G 
df* yield \g{u) - giU.)\ < ||R.,*|U*||m - Ufl\x. ’ ■ 

1.3. Outline. In Section 2, we propose three algorithms which are analyzed in the 
sections below, define the used mesh-refinement strategy, and outline the main result. 
Moreover, we provide the abstract framework in terms of four axioms for the estimators. 
Section 3 proves optimal convergence rates for each adaptive algorithm. In Section 4 
we apply the abstract theory to conforming goal-oriented FEM for second-order elliptic 
PDEs. Section 5 covers goal-oriented EEM for the evaluation of some weighted boundary 
flux, whereas Section 6 applies the abstract theory to goal-oriented adaptivity for BEM. 
The final Section 7 discusses our results and points at extensions and open questions. 

2. Adaptive Algorithms for the Estimator Product 
W e consider an adaptive algorithm which allows to drive the estimator product 

(8) Vu,*Vz,* ■= {'^ VuATf) ' ( XI ' 

rer* rer* 

to zero with optimal rate. This includes, in particular, the problem class from Section 1.2, 
but also covers adaptive BEM for the approximation of point errors; see the recent own 
work [20]. We suppose that each admissible triangulation (see Section 2.2 below) allows 
for the computation of the error estimators rju,* and where the local contributions 
are (at least heuristically) linked to the elements T G 7^, cf. (8). To abbreviate notation, 
we shall write 

Vw,* ■= VwAr), w,*(^*) •= ( X ^ for tc G {u, z} and all Z4 C %. 

Teu„ 

2.1. Adaptive algorithm. We consider two adaptive algorithms, which have been 
proposed and analyzed in [37] (Algorithm A) and [7] (Algorithm C) for goal-oriented 
adaptive FEM for the Poisson problem, and propose a slight modification of the algorithm 
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from [37] (Algorihm B) which is also related to that of [29]. Note that all algorithms 
differ only in the marking strategy: Algorithms A-B employ a separate Dorfler marking 
in step (iii)-(iv), whereas Algorithm C employs a combined Ddrfler marking in step (iv). 

Algorithm A. Input: Initial triangulation To, marking parameter 0 < 6* < 1, and 

C*mark ^ 1 • 

Loop: For all £ = 0,1,2,3,... do (i)-(vi).- 

(i) Compute refinement indicators rju/iT) for all T G Tr- 

(ii) Compute refinement indieators rjz^fiT) for all T ETe- 

(hi) Determine a set Aiu,e Q Te of up to the multiplicative factor Cmark minimal car¬ 
dinality such that 

(9) 0 < riu,fiMu,if ■ 

(iv) Determine a set A4z,e QTi of up to the multiplicative factor Cm ar k minimal car¬ 
dinality such that 

(10) Or)l,><r]zA^z,ef- 

(v) Choose Aii G {Mu,£, fidzA to be the set of minimal cardinality. 

(vi) Let Ti+i := ref ine(7r, Adr) be the coarsest refinement of% such that all marked 
elements T G Air have been refined. 

Output: Sequence of successively refined triangulations % and corresponding error es¬ 
timators r]u,£, Vz,£ for all £ E No. ■ 

Remark 2. In the frame of Section 1.2, the computation ofrju/ in (i) and rjz/ in (ii) usu¬ 
ally requires to solve the primal (2) and the dual problem (5) to obtain the approximations 
Ui resp. Z(. ■ 

Remark 3. To construct a set A4„/ with minimal eardinality (i.e., Cm ar k = 1) which 
satisfies the Ddrfler marking (9), one needs to sort the refinement indieators. This results 
in logarithmie-linear costs. For Cmar k = 2, the algorithmie eonstruetion of some set A4.u,£ 
which satisfies (9), is possible in linear complexity [42], ■ 

Next, we propose a modified version of Algorithm A which allows for more aggressive 
marking in step (v), i.e., less adaptive steps. A similar but non-optimal algorithm has 
been proposed and tested in [29]. 

Algorithm B. Input: Initial triangulation To, marking parameter 0 < 0 < 1, and 
Cmark; > 1. 

Loop: For all £ = 0,1,2,3,... do (i)-(vi).' 

(i) Compute refinement indicators rjuAT) for all T G 7r- 

(ii) Compute refinement indicators rjz/{T) for all T G 7r- 

(iii) Determine a set Aiu,£ Q T of up to the multiplicative factor Cmar k minimal car¬ 
dinality such that 

( 11 ) 0 < rjuA^u/f ■ 

(iv) Determine a set A4z,£ of up to the multiplicative factor Cmark minimal car¬ 
dinality such that 

( 12 ) Ogl^<gzA-^z,£f- 

(v) Choose G {A4u/, to he the set of minimal eardinality and ehoose Adf C 

^ U J^z/ SUeh that ^ 
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(vi) Let %+! := ref ine(7£, be the coarsest refinement of% such that all marked 
elements T G Aii have been refined. 

Output: Sequence of successively refined triangulations % and corresponding error es¬ 
timators riu/, riz,i for all i £Nq. ■ 

Remark 4. In our numerical experiments below, we choose M.i as follows: Having picked 
Aie to be the minimal set amongst A4.u,i ond M.z,e, we enlarge Aie by adding the largest 
fi^Aii elements of the other set, e.g., if ffAiu,i < then A4( consists of AAu,i plus 

the ifAiu,e largest contributions of AAz/- This yields ~ 2. ■ 

Remark 5. In [29], the authors consider Algorithm B, but define AAi := AAu/ U AAz,£ 
in step (v). While this also leads to linear convergence in the sense of Theorem 16, 
[29] only proves suboptimal convergence rates min{s, t} instead of the optimal rate s -\-t 
in Theorem 20; see [29, Section 4], We note that the strategy of [29] leads to linear 
convergence riu,e+n < Cq^riu/ und riz,i+n < Cq’^Pz/ for either estimator and all i,n E No, 
where (7 > 0 and 0 < q < 1 are independent constants, while the optimal strategies 
considered in this work only enforce linear convergence rju,e+nOz,e+n < Cq^rju/Oz/ for the 
estimator product. ■ 

Finally, the following algorithm has been proposed in [7]. 

Algorithm C. Input: Initial triangulation To, marking parameter 0 < 6 < 1, and 
Omark ^ 1 • 

Loop: For all i = 0,1, 2,3,.. . do (i)-(vi); 

(i) Compute indicators rju,e{T) for all T eTi. 

(ii) Compute indicators rjz/iT) for all T E %. 

(hi) Assemble refinement indicators Pe(Tfi := ^ + rfi^ fr}z,e.{TY for all T ETe. 

(iv) Determine a set AAi F % of up to the multiplicative factor Omark minimal cardi¬ 
nality such that 

(13) OpHpeiMef. 

(v) Let 7r+i := ref ine(7r, Air) be the coarsest refinement ofTe such that all marked 
elements T E AAi have been refined. 

Output: Sequence of successively refined triangulations Te and corresponding error es¬ 
timators r]u,i, r}z,i for all £ eNq. ■ 

Remark 6. Building on the results o/[37], it is claimed and empirically investigated in [7] 
that the combined Dorfler marking (13) requires less adaptive steps to reach a prescribed 
accuracy. We note that this has only been proved rigorously in [7] for the Poisson problem 
with polynomial data, while general data have to satisfy a certain saturation assumption 
for the related data oscillation terms; see [7, eq. (4.4)] and [7, Theorem 4.1], We note 
that this assumption also restricts the quasi-optimality analysis of [7] which, in the spirit 
of [37], relies on the (constrained) contraction of the total error. ■ 

2.2. Mesh-refinement. We suppose that the mesh-refinement is a deterministic and 
fixed strategy, e.g., newest vertex bisection [43]. Unlike [7, 37], we do not require the 
interior node property guaranteed by bisec5-refinement of marked elements [39]. For 
each triangulation T and marked elements A4 C 7~, we let T := refine(T, Af) be 
the coarsest triangulation, where all elements T E A4 have been refined. This may, 
in particular, include the preservation of conformity or bounded shape regularity. We 
write T E refine(T), if there exist finitely many triangulations T^^\ ... and sets 
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>f(i) C such that T = V = and = refine(r(^-i),for 

all j = l,...,n, where we formally allow n = 0, i.e., T = 7”^°^ G refine(7”). To 
abbreviate notation, let T := refine (To), where To is the given initial triangulation of 
Algorithms A-C. 

2.3. Main result. Our main result requires the following abstract approximation 
class for the error estimator r]u,e resp. Let T^v := {T G T : #T— #7o < A"} denote 
the (finite) set of all refinements of To which have at most N elements more than To- For 
s > 0 and w G {«, z}, we write w G if 



where is the error estimator associated with the optimal triangulation % G T^v- In 
explicit terms, ||w||as < oo means that an algebraic convergence rate 0{N~^) for the 
error estimator is possible, if the optimal triangulations are chosen. 

For either algorithm, our abstract main result is twofold: First, we prove linear con¬ 
vergence (Section 3.1): For each 0 < g < 1 there exists some n such that for all £ G N, 
the reduction of riu/ by the factor q requires at most n steps of the adaptive loop, i.e., 
Vu,£+nVz,e+n < QVu/Vz/- Secoud, we prove optimal convergence behavior (Section 3.3- 
3.5): With respect to the number of elements N ~ #Tr —#7o, the product rju,iTjz,£ decays 
with order £T(A“(^+*)) for each possible algebraic rate s + t > 0, i.e., ||u||as + IklU* < oo. 
This means that the adaptive algorithms will asymptotically realize each possible alge¬ 
braic decay. 

Remark 7. Since our analysis works with the estimator instead of the error, it avoids 
the use of any (discrete) efficiency bound. Compared to [7, 37], this allows to use simple 
newest vertex bisection instead of bisecS-refinement for marked elements. Lemma 22 below 
states that for standard FEM our approximation class coincides with that of [7, 18, 37] 
which is defined through the so-called total error (i.e., error plus data oscillations). ■ 

2.4. Axioms of Adaptivity. In the following, let dl^„(-, •) : T x T ^ M>o denote a 
distance function on the set of admissible triangulations which satisfies 


, r") < di,,(r, T) + d^{r , r") for aii r, r, r" g t, 

d^{T,r) < C'distdl^(r, T) for all T, T' G T, 


with some uniform constant Cdist > 0; see also Remark 11 below. 

The convergence and optimality analysis of Algorithm A is done in the frame of the fol¬ 
lowing four axioms of adaptivity [15], where axiom (A3) is slightly relaxed when compared 


to [15]: 


(Al) Stability on non-refined elements: There exists a constant (Tstb > 0 such that for 
each triangulation Tr G T and all refinements 7) G refine(Tr) the corresponding 
error estimators satisfy 


lw,*(TrnT;) (7^ ft;)I < C'stbdl«,(Tr ,t;) 


(A2) Reduction on refined elements: There exist constants 0 < gred < 1 and Cred > 0 
such that for each triangulation Tr G T and all refinements Tt G ref ine(Tr) the 
corresponding error estimators satisfy 


VwffiTffiTif < gred VwMXKf + C'red d^{T£ , 
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(A3) Quasi-orthogonality: Let Ti^ be the (possibly finite) subsequence of triangulations 
Ti generated by Algorithm A, B, or C which satisfy the Dorfier marking on the 
refined elements, i.e., 

(14) < Vw/niTeATin+if- 

Then, for all s 0, there exists some constant C^Qj’h}i(<s^) 0 such that for all 

n < for which 7^^,..., are well-defined, it holds 

N 

(dl^(7^,+i < CoM£)vY- 

j=n 

(A4) Discrete reliability. There exists a constant Crei > 0 such that for each triangula¬ 
tion 7^ G T and all refinements 71^ G refine(7^), it holds 

d«)(7X,7r) < Cie\'r)w,e{'^w{T£,%)), 

where the set TZw{%,%) satisfies %\7l C TZw{%,%) ^ % and #7^«;(7r,7X) < 
C'rei #(7r\71^) and hence consists essentially of the refined elements only. 

We recall some elementary observations of [15]. 

Lemma 8 (quasi-monotonicity of estimator [15, Lemma 3.5]). There exists a constant 
Cmou > 0 which depends only on stability (Al), reduction (A2), and discrete relia¬ 
bility (A4), such that for all % E T and all refinements 71 G refine(7^), it holds 
Vw,* — ^yi^onVw,£' ® 

Lemma 9 (optimality of Dorfier marking [15, Proposition 4.12]). Suppose stability (Al) 
and discrete reliability (A4). For all 0 < 6 < 6^, := (1 -|- there exists some 

0 < «:* < 1 such that for all Te ET and all refinements 71 G ref ine(7£), it holds 

( 15 ) vl,*< 

where 7Zw{Te.,%) is the set of refined elements from (A4). ■ 

Remark 10 (validity of quasi-orthogonality for £ = 0). Suppose that (A1)-(A4) are valid. 
Arguing along the lines o/[15. Proposition 4.11], one then sees that quasi-orthogonality (A3) 
holds even with s = 0 and (Porth depends only on (Pconv? Qconv from Proposition 13 and 
Crel- 

Remark 11. (i) In the setting of Section 1.2, let w G {u,z} with IT* G {T*, Z*} being 
the corresponding Galerkin solution for % eT. The abstract distance is then defined by 
, Z) := a(lT* - IT,, IT* - IT,)^^ ~ ||1T* - 1T,|U; see Section f-6 below. 

(ii) Suppose that the bilinear form «(•,•) is additionally symmetric, and let |||n||| := 
a(n,n)^/^ denote the induced energy norm which is an equivalent norm on X. Then, 
nestedness C C X^ of the discrete spaces for all k > m > n implies the Galerkin 
orthogonality 

lIlTfc - VTnf + IliVTn - iwf = lIlTfc - iwf for allk>m>n. 

This and (A4) imply 

N N 

j=n j=n 

< < vl,,. 

and hence the quasi-orthogonality (A3) with £ = 0. ■ 
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2.5. Generalized linear convergence. The following estimator reduction is first 
found in [18] and, e.g., proved along the lines of [15, Lemma 4.7]. Since we need a slightly 
stronger result than that of [15], which covers arbitrary refinements G refine(7£+i) 
instead of Ti+\ only, we include the proof for convenience of the reader. 

Lemma 12 (generalized estimator reduction). Let 0 < d < 1. Let 7^ G T and Te+i G 
refine(7^) be given triangulations and suppose that the refined elements satisfy the 
Dorfler marking 

( 16 ) Oril^^<ri^^i{Te\%+if. 

Then, there exist eonstants 0 < gest < 1 Cest > 0 whieh depend only on assump¬ 
tion (A1)-(A2) and 9 such that for all % G ref ine(7^+i), it holds 

(17) vl,* < Qest vl,e + C'est d[^„(7;, Tef. 

Proof. For each 5 > 0, the Young inequality (a + bfi < (1 + 6 ) 0 “^ + (1 + S~^)b‘^ and 
stability (Al) yield 

n TeY < (1 + (^) 9w,iO~k n TiY + + (^ ^) dyj{7l , TiY 

= (1 + (5) ryi, - (1 + 5) Vu,Ari\%f + + r') d^(r* , Tef. 

Together with reduction (A2), we obtain 

7^* = VwA'^* 61 Tif + 9wA'T*\Tif 

— (1 + Vw,e ~ (1 + ~ gred) 9wA'T~Aj~kY + (61stb(l + <^ ^) + C'red) dltt,(71( , T~ff. 

With %\%: 5 TifTi+i and the Dorfler marking (16), we see 

VwA'^A'T*^ > VwA'^A%+if > 

Combining the last two estimates, we prove (17) with Cest = C'stb(l + + ^^red and 

gest = (1 + (5) — (1 + (^ — gred)^ ^ 1 — (1 — gred)^ < 1 as (5 —?• 0. For Sufficiently small h > 0, 
it thus holds 0 < gest <1- B 

The following result generalizes [15, Proposition 4.10] to the present setting. We note 
that (A4) enters only through the quasi-monotonicity of the estimator (Lemma 8). 

Proposition 13 (generalized linear convergence). Let % be a sequence of successively 
refined triangulations, i.e., Te G ref ine(7£-i) for all £ G N. Let 0 < 9 < 1. Then, there 
are constants 0 < gconv < 1 and Cconv > 0 which depend only on (Al)-(A4) and 9, such 
that the following holds: Let i,n & Nq and suppose that there are at least k < n indices 
i<ii<i 2 <-- - < A such that the Dorfer marking (16) is satisfied on the 

refined elements, i.e., 

(18) 9 r]lj^ < (Te^ \Te^+iA for all j = I,... k. 

Then, the error estimator satisfies 

(19) 9w,t+n — Aconv Qconv Vw,^- 

Proof. To abbreviate notation, set jo := i. Note that G ref ine(7£j.+i). Therefore, 
the estimator reduction (17) shows for all £ > 0 and all 0 < g < 

k k 

Y. ^ I] ( QestVw/i ^ 

i=k—j i=k-j 

k 

~ ^ ^ + ^est^)Vw/i 5 7ii) ~ ’ 


Choose £ < (1 — q'est)C'ggt so that k := 1 — (gest + Cest^) > 0. Use (A3) to obtain for all 
0<j<k 


( 20 ) 


^ ^ ^ Vw/i^i — '^w/k-j ^ V ^Tii) 

i=k—j i=k—j 

^2 


< (1 + C'estC'orth(£))?7t„,4 . 



With (7 := (1 + C'estC'orth(£))/«^ > C mathematical induction below shows 


k 

(21) vl,e, < (1 - 0 < i < fc. 

i=k-j 

To see (21), note that the case j = 0 holds with equality. Suppose that (21) holds for 
j < k. This induction hypothesis and (20) show 


'nl 


,4 < (1 - C y = (1 - C' ')^'(( J] <q) - 7i4-0+l)) 

i=k-j i=k-{j+l) 

(20) JL^ 

<(i_c-4‘ 4,„ 

i=k-{j+i) 


which proves validity of the induction step. Hence, the assertion (21) holds for all j < k. 
By use of (21) for j = fc — 1 and (20) for j = k, we obtain 


( 21 ) 


CmlnVl,i+n < vIa < (1 “ C' ^ < (1 “ C' 


2 


( 20 ) 


1=1 

y-l\k-lr^ ^2 


1=0 


< (1 - c-r-^CvU = (1 - c-^c/ii - C-^) vIa 

This concludes the proof with Cconv = C'C'mon/(l — C~0 and gconv = (1 “ C~0- 


3. Optimal Convergence of Adaptive Algorithms 

Throughout this section, we suppose that the error estimators riu,i and riz/ satisfy the 
respective assumptions (A1)-(A4) of Section 2.4. Without loss of generality, we suppose 
that riu,e and rjz/ satisfy the axioms (A1)-(A4) with the same constants. 

Remark 14. The axioms (A1)-(A4) are designed to eover weighted-residual error es¬ 
timators in the frame of FEM and BEM. However, as is shown in [15, Section 8] for 
optimal adaptivity for the energy error, it is sufficient that for w e {u, z} the error es¬ 
timator r}w,e used in the adaptive algorithm is locally equivalent to some error estimator 
rjyj^i which satisfies (Al)-(A4), i.e., 

^ AwMA) AUT) ^ for all T e Ti, 

where u)(,{T) denotes a patch ofT. Then, the convergence (Theorem 16) as well as opti¬ 
mality results (Theorem 11, 20, and 21) remain valid. lUe leave the details to the reader, 
but note that such arguments cover averaging-based error estimators, hierarchieal error 
estimators, as well as estimators based on equilibrated fluxes; see [15, 33]. ■ 


3.1. Linear convergence of Algorithms A, B, and C. The following result gives 
a criterion for linear convergence which is satisfied for either adaptive algorithm (Algo¬ 
rithm A, B, and C). As a consequence, linear convergence is independent of Cmark, and 
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we may formally also choose Cmark = oo. Finally, linear convergence (22) does only rely 
on discrete reliability (A4) to ensure quasi-monotonicity of the estimator (Lemma 8). In 
the frame of the Lax-Milgram lemma from Section 1.2, the quasi-monotonicity already 
follows from classical reliability (6); see [15, Lemma 3.6]. 

Proposition 15. Let the sequenee of meshes (7^)^eNo C T satisfy Ti+i G refine(7^) for 
all £ G Nq. For all i G No, let the set F4e := 7^ \ 7e+i satisfy either (9) or (10). For all 
0 < 0 < 1, there are constants 0 < giin < 1 and C\i^ > 0 which depend only on (Al)-(A4) 
and 9, such that there holds linear convergence in the sense of 

(22) Vu,i+nVz,(+n — C*iin5iiii%,rhz,r fof' 0,11 £,n ^ No- 

Proof. The assumptions on (7 ^)£gNo imply that for n successive meshes Tj, j = i,, i+n, 
Tj\Tj+i satisfies fc-times the Dorfier marking for riu/ and (n —fc)-times the Dorfier marking 
for riz/. According to Proposition 15, this implies 

vl,e+n < C'conv gconv vl,e ^S Well aS < Cconv Cnv vlf 

Altogether, this proves 

2 2 ^ k 2 2 

Vu,(+n '9z,i+n — ^conv 7conv 9u,i 9z,(' 

This concludes (22) with gun = qlluv and = (Pconv ® 

Theorem 16. For all t) < 6 < 1, there are constants 0 < giin < 1 and Ciin > 0 which 
depend only on (A1)-(A4) and 9 such that Algorithms A, B, and C are linear convergent 
in the sense of (22) 

Proof for Algorithm A. In each step of Algorithm A, the set satisfies either the Dorfier 
marking (9) for r]u,e or (10) for pz^i- Since Adi C the assumptions of Proposi¬ 
tion 15 are satisfied. This concludes the proof for Algorithm A. ■ 

Proof for Algorithm B. Analogously, the set Adi satisfies in each step of Algorithm B 
either (11) or (12). Since Ad(, C A4^ Q Ti \ Te+i, Proposition 15 applies and concludes 
the proof for Algorithm B. ■ 

Proof for Algorithm C. Note that = 2 r]^^r 7 ^^. Therefore, (13) becomes 

In particular, this shows that 

0 vl,e < VuA^ef or 9 < rizAMef ■ 

Since Adi ^% \ 7r+i, Proposition 15 applies and concludes the proof. ■ 

3.2. Fine properties of mesh-refinement. The following Theorem 17 states optimal 
convergence behavior. Unlike linear convergence, the proof of optimal convergence rates 
is more strongly tailored to the mesh-refinement used. First, we suppose that each refined 
element has at least two sons, i.e., 

(23) #(T\r') + #r < AT' for all r G T and T G ref ine(r). 

Second, we require the mesh-closure estimate 

(24) AT - AT ^ enmesh E AAdj for all £ G N, 

i=o 

where Cmesh > 0 depends only on To. This has first been proved for 2D newest vertex 
bisection in [11] and has later been generalized to arbitrary dimension d > 2 in [43]. 
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While both works require an additional admissibility assumption on To, this has at least 
been proved unnecessary for 2D in [32], Finally, it has been proved in [18, 42] that newest 
vertex bisection ensures the overlay estimate, i.e., for all triangulations T, T' G T there 
exists a common refinement T ®T' £ ref ine(7~) fl ref ine(7~') which satisfies 

( 25 ) #(r © T) <ihT +#r - #ro. 

Although not used explicitly, we note that for newest vertex bisection, the triangulation 
T (BT' is, in fact, the overlay of T and T'. For ID bisection (e.g., for 2D BEM compu¬ 
tations in Section 6), the algorithm from [2] satisfies (23)-(25) and guarantees that the 
local mesh-ratio is uniformly bounded. For meshes with first-order hanging nodes, (23)- 

(25) are analyzed in [12], while T-spline meshes for isogeometric analysis are considered 
in [38]. 

3.3. Optimal convergence rates for Algorithm A. Our proofs of the following 
theorems (Theorem 17, 20, 21) follow the ideas of [37] as worked out in [20]. We include 
it here for the sake of completeness and a self-contained presentation. 

Theorem 17. Suppose that the mesh-refinement satisfies (23) as well as the mesh-elosure 
estimate (24) and the overlay estimate (25). Let 0 < 0 := (1 + he 

sufficiently small. Then, there exists a constant Copt > 0 which depends only on 0, 
Omesh, o-nd (Al)-(A4), such that for all > 0 the assumption {u,z) G x Aj implies 
for all £ G No 

(26) - g ,\ < Vu . eri..e < II“IU.II^IIa. 

i.e., Algorithm A guarantees that the estimator product decays asymptotically with any 
possible algebraic rate. 

Corollary 18. Assume that the estimators both have finite optimal eonvergenee rate, i.e., 
-Smax := sup{s > 0 : ||u||as < Oo} < OO and tmax := sup{t > 0 : ll^llAt < 00} < 00. 
Then, (26) implies for all 0 < s < s^ax o,nd 0 < t < tmax and a// £ G N 
Wr < (#7^ - «)-* as well as 

where the hidden constants additionally depend on s^ax — s > 0 resp. t^ax — t > 0. ■ 

The heart of the proof of Theorem 17 is the following lemma. 

Lemma 19. For any 0 < d < := (1 + and i G No, there exist eonstants 

Cl, (7*2 > 0 and some refinement % G refine(77) sueh that the sets 77„(77,77) cind 
from the discrete reliability (A4) satisfy for all s,t>0 with {u,z) G A^ x Aj 

(27) max{#7^„(77,r*), < C, (C'2||u||Ajkl|Aj'/^*+*^ 

Moreover, either or 7lz{Te,7l) satisfies the Ddrfler marking, i.e., it holds 

(28) < r]u,i{Tiu{%,%)Y or < r]z,i{F-z{Ti,%)Y. 

The constants Ci,C 2 depend only on 9, (Al),(A2), and (A4). 

Proof. Adopt the notation of Lemma 9. For e := the quasi-monotonicity 

of the estimators (Lemma 8) yields £ < < IlifllAjl^^llAt < oo- Choose the 

minimal A G No such that ||M||Aj|'2(||At < £ (A + l)®+h Then, Lemma 8, the definition of 
the approximation classes, and the choice of A give 

9u,*9z,-k — ClY^onVu,elVz,e2 — Cinon(A^ © 1) ^ ^ ||^||As ||^||Af © Cmon^ /^*Vu,£Vz,£- 
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This implies or 77^* < and Lemma 9 hence proves (28). It remains 

to derive (27). First, note that 

(29) max{#77„(7^, T*), #77,(7^, T*)} C^ei #(7^\r*) 0^#% - #7^), 

since refined elements are refined into at least two sons (23). Second, minimality of N 
yields 

with C := (||M|UJk||At)^/^*'^*nC'moii'«*)“^^^*'^*^ = (C'mon/«7Ml“IUJk||At)^/^*'^*^- ChoOSC 
7^1,7^2 G Tjv with r7„,£j = minr^eT^v Vu,* and 77^,= minr^eTjv Vz,*- Define 
and 71 := Te (S)Te. The overlay estimate (25) yields 

( 25 ) ( 25 ) 

(30) #r* -#% < #% - #ro < + ihTe, -2#%<2N< 2C (77^,,77,,,)-'/(*+*). 

Combining (29)-(30), we conclude (27) with Ci = 2(7161 and C 2 = Craon/t^*- B 

Proof of Theorem 1 7. According to (28) of Lemma 19 and the marking strategy in Algo¬ 
rithm A, it holds for all j G No 

(31) = YCill\.\^AA.^j , TT^Adxj} ^ Cmark 12iax-[77^:77.u(7j, 7*) , fj^'PzifPj ■,7f)'\ ■ 

With the mesh-closure estimate (24) and estimate (27) of Lemma 19, we obtain 

( 24 ) ( 27 ) 

#7^-#ro < C'^eshJ]#Ad, < C^eshC„,arkCi(C2||u||Ajk||Aj'/(^+*) 

7=0 7=0 

Linear convergence (22) implies 

< C'lin gfin^77^,77^,7 for all 0 < 7 < £ 

and hence 

With 0 < g := < 1, the geometric series applies and yields 




i-i 


j=0 j=0 

Combining this with the first estimate, we obtain 


< 


r 


1/(s+t) 


lin 


1_ 

^lin 


{'nu,er]^ 


z,£ 


l/(s+t) 


#7^ - #ro < 


enmesh C^mark C i 
1_ 

^lin 


(ainC2 ||M||Ajkl|Aj'/(^+*) 


Rearranging this estimate, we conclude (26) with (7opt = max{Ciin(72, (7mpsh( 7mar k(7i}. 


3.4. Optimal convergence rates for Algorithm B. 

Theorem 20. Let d* := (1 + (7stbC'rei)~^. For any 0 < d < 6**, Algorithm B guarantees 
optimal algebraic convergence rates in the sense of Theorem 17 and Corollary 18. 

Proof. We note that Lemma 19 is not affected by the marking strategy and hence remains 
valid. To conclude the proof, we only need to show that estimate (31) also remains true. 
Since A4j in Algorithm B is a set of minimal cardinality up to the factor (7()^ark^mark 
which satisfies either (9) or (10), estimate (31) holds with different constants, i.e., 

ifMj < (7(„arkC'markmax{#7^„(7),7;), #n^{Tj,%,)}. 

Therefore, the claim follows with (7opt = max{(7ii„(72, (7meshC'(jjark^markC'i}. ■ 
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3.5. Optimal convergence rates for Algorithm C. 

Theorem 21. Let := (1 + C'stbC'rei)”^- For any 0 < 9 < ^*/2, Algorithm C guarantees 
optimal algebraic eonvergenee rates in the sense of Theorem 17 and Corollary 18. 

Proof. We only need to show that (31) remains valid. Note that 0 < 26* < 6**. Therefore, 
estimate (28) of Lemma 19 yields 


< r]u,j{F.u{Tj,Z)y or 207]^ < . 



According to the marking strategy in Algorithm C, we obtain 


which is (31). Therefore, the claim follows with Copt = niax{C'iinC' 2 , C'meshC'markC'i}- ® 

4. Goal-Oriented Adaptive FEM for Second-Order Linear Elliptic PDEs 

In this section, we extend the ideas from [23] and prove that our abstract frame of con¬ 
vergence and optimality of goal-oriented AFEM applies, in particular, to general second- 
order linear elliptic PDEs. 

4.1. Model problem. On the bounded Lipschitz domain D C and for given 
fi,gi G L?‘{Q) and / 2,^2 ^ L‘^{ny, we aim to compute 



Jn 

where u G Hl{Q) is the weak solution to 

(32) Cu := —div(AVM) -|- b • Vu -|- cm = /i -|- div /2 in fl and u = 0 on P := dfl. 

For all X E Ll, A{x) G Mgym ^ symmetric matrix with A G 1F^’°°(D; Moreover, 

b{x) G is a vector with b G 1F^’°°(D; M*^), and c{x) G M is a scalar with c G L°°(D). To 
formulate the residual error estimators in (35)-(36) below, we additionally require that 
div/ 2 ,divgr 2 exist in elementwise on the initial mesh To and that the edge jumps 

satisfy [/2 • n], [^2 ' ^ L?‘{dT) for all T E %. (These assumptions are for instance 

satisfied if are 7o-piecewise constant.) Note that C is non-symmetric as 

(33) Cw 7 ^ = —divAVtc — b • Vin + (c — divb)w. 

We suppose that the induced bilinear form 

a{u, v) := {Cu ,v)=[ AVu • Vn + 6 • Vuv + cuv dx for u,v E X := H^{n) 

Jn 

is continuous and iTg (r2)-elhptic and hence fits in the frame of Section 1 . 2 . The right-hand 
side of ( 1 ) reads f{v) := fiv — f 2 • Vv dx. 

4.2. Discretization. For a given regular triangulation of D and a polynomial degree 
p > 1, define := {V E LF‘{Vt) : V\t is polynomial of degree < p for all T E %}. 

We consider A* := Sfilf) \= V^ilf) fl Hl{Q) and let G A* be the unique FEM 

solutions of (2) resp. (5), i.e., 


(34a) U,eSI{%) such that a{U,,V,) = f{V,) for all F* G 5^(7;), 

(34b) Z,eSI{%) such that a(W, Z*) = p(W) for all W G 5^(7;). 
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4.3. Residual error estimator. For T ETk, let hx '■= and C\tV := —div|TA(W)+ 

b • W + cV. Then, the residual error-estimator for the discrete primal problem (2) reads 

(35) := hlWClrU. - - div+ AtIIKAVC/. + f,) ■ n]\\%f„^y 

The residual error-estimator for the discrete dual problem (5) reads 

(36) ?7^,*(T)2 := h%\\L^\TZ^ - gi - div^2lli2(r) + ^t||[(AVZ* + ^ 2 ) ' 

where C^\tV := -div|TA(Vl^) - 6 ■ VF + (c - div6)l^. 

The error estimators satisfy reliability (6); see, e.g., [1, 44], The abstract analysis of 
Section 1.2 thus results in 

(37) \g{u) - g{U^)\ < 

and we aim for optimal convergence of the right-hand side. Moreover, efficiency and 
the Cea lemma prove that the estimator based approximation class A* from Section 2.3 
coincides with the approximation class based on the total error used, e.g., in [7, 18, 37]. 

The following result is proved in [23, Lemma 5.1] for /2 — 0 = g 2 , but holds verbatim in 
the present case. 

Lemma 22. Let w G {u, z}. There holds w £ Ag if and only if 

sup ((N + iy min ( min ||w - V7|U + osc^,*(14)) ) < (X), 

A/'gNo T^GTiv / 

where osc^„,*(l/*)^ = I^Ter* Iff and 

osclJT, K.) - - Yif-^CWV, - - div/2)||i.,„ 

+ AtII (1 - n%-')l{AVK + h) ■ n]|li.,»Tn!!). 

osc=_.(r, K.) - 411(1 - n|>’-")(£^|rK, - Ji - divg^dli.ij,) 

+ AtIKI - n%-')l{AVV, + g,) ■ nllli.iarnn)- 

Here, Wp : Lf{T) V‘^{T) denotes the L‘^-orthogonal projection onto polynomials of 
degree q and 11^^ : L‘^{dT) ‘P'^{Sqt) denotes the -orthogonal projeetion onto (dis- 

eontinuous) piecewise polynomials of degree q on the faces ofT. ■ 

4.4. Verification of axioms. With newest vertex bisection from [43] as mesh- 

rehnement strategy, the assumptions of Section 3.2 are satished. It remains to verify 
the axioms (Al)-(A4), where d^(7^ , %) := a{We - W*, W, - ~ ||1W - lW||^i(a) 

and Wi resp. IF* are the corresponding FEM approximations of w G {u, z}. 

Theorem 23. Consider the model problem of Section f.l. Then, the conforming dis¬ 
cretization (34) with the residual error estimators (35)-(36) satisfies stability (Al), re¬ 
duction (A2) with Qred = 2“^/'^, quasi-Orthogonality (A3), and discrete reliability (A4) 
with TZyy{T£,%) = %\T-k and w G {u,z}. In particular, the Algorithms A-C are linearly 
convergent with optimal rates in the sense of Theorem 16, 17, 20, and 21 for the upper 
bound in (37). 

Proof of Theorem 23, (A1)-(A2) and (A4). The work [18] considers some symmetric model 
problem with 6 = 0 and c > 0 as well as /2 = 0 = ^ 2 - Stability (Al) and reduction (A2) 
are essentially part of the proof of [18, Corollary 3.4]. The discrete reliability (A4) is 
found in [18, Lemma 3.6]. Both proofs transfer verbatim to the present situation with a 
non-symmetric differential operator and general / 2 , 3^2 • ® 
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Lemma 24. In the setting of Theorem 23 and for Algorithm A-C, the Galerkin approx¬ 
imations Ui and Zi converge in the sense of 

(38) lim \\Uoo - Ue\\m{n) = 0 = lim \\Zoo - Ze\\m{n), 

for certain Uoo, Z^o G Moreover, there holds at least Uoo = u or Z^o = z. 

Proof. Note that adaptive mesh-refinement guarantees nestedness C A”* for all 7^ G T 
and %: G refine(7r). As in [15, Section 3.6] or [5, Lemma 6.1], the Cea lemma thus 
implies a priori convergence in the sense that there exist Uoo,Zoo G AA, := u 


IIq(Q) such that 


'^eNo — 


lim \\Uoo — — 0 — lim ||Zoo — 


For w G {u,z}, let denote the subsequences which satisfy 

for all n G M. 

There holds if{Iw,n '• n G N} = (X) for at least one w G {u,z}. While this is obvious for 
Algorithm A and Algorithm B, it follows for Algorithm C from the proof of Theorem 16. 
For this particular w, the estimator reduction from Lemma 12 reads 

+ C'estdl«,(7^,„.n+1 ,for all n G N. 

The a priori convergence (38) implies ^ 0 as n —>■ oo. Elementary 

calculus thus yields lim^^oo hw.c.n = Oj [13, Corollary 4.8] resp. [4, Lemma 2.3]. 

Reliability (6) of 77^^ concludes lim^^oo = 0, i.e. w = Woo- ■ 

Proof of Theorem 23, (A3). Recall the sequence from (A3). With the a priori conver¬ 
gence of Lemma 24, the proof of [23, Lemma 3.5] applies and shows the weak convergence 
in H^(Q) for Woo G {Uoo, Zoo} 


Ifooo - 


0 and 




n+1 


Wa 


0 as £ —>■ 00. 


\\Wo,-wam») ” inu„ - iuJIhmq) 

Define d^wiToo, •) := a(Woo — (O^Hoo — With this, [23, Proposition 3.6] applies 

for the primal as well as the dual problem and shows that given any h > 0, there exists 
js G N such that all j > js satisfy 

,Tl,f < jf-rdwiXo ,Tt,f - 


(39) 


The discrete reliability (A4) and the convergence (38) yield 
(40) dl^„(7;o ,7^,-) = lim dl^(7^j, < C,e\r)w,ey 


k^oo 


With (39)-(40), the quasi-monotonicity from Lemma 8 (since (Al), (A2), (A4) have 
already been verified) implies for h = 1 — 1/(1 -|- sC~^) and hence 1/(1 — h) = 1 -|- eC^ 


1-2 

rel 


that 


N 


Yi (d„(y„ ,T,,f - EC-^d4T„ .T,p) 


j=n 


( 39 ) 

- “ ^C;2)d^{roo , Ti^f - d^ifToo , 


(41) 


J=J5 


j=n 


jd-l 


( 40 ) 


< d^('7^ , )^ + C/gi ^ < (1 js)C^eiC, 




j=n 
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Figure 1. Example from Section 4.5: The initial mesh To (left) and the 
triangles Tf (bottom left) and Tg (top right) indicated in gray. An ap¬ 
proximation to the primal solution (middle) and dual solution (right) on a 
uniform mesh with 256 elements, where the singularities of both are clearly 
visible. 


Another application of the reliability (40) shows 

^ ( 40 ) ^ 

j=n j=n 

( 41 ) 

— (It 2<5)^rel^monTiy/„- 

This proves (A3) with C'orth(£) := (I + js)C^eiCmon- ■ 


4.5. Numerical experiment I: Goal oriented FEM for the Poisson equation. 

We consider a numerical example proposed in [37, Example 7.3] for the Laplace operator 
in 2D, while a nonsymmetric second-order elliptic operator is considered in Section 5.5. 
The goal of this first experiment is to verify the optimal convergence of Algorithm A- 
C as predicted by theory, and to compare the various algorithms as well as standard 
AFEM (i.e., non-goal-oriented adaptive FEM, where Aii A4u,e resp. A4^ := in 
Algorithm A); see, e.g., [15, 18, 23, 42]). 

We consider the Poisson model problem (i.e., A — I, b = 0, and c = 0) on the 
unit cube fl = (0,1)^ C Unlike [37] which considers quadratic elements p = 2, 
we apply cubic elements p = 3 (unless stated otherwise). The initial mesh To is shown 
in Eigure 1 (left), where also the triangles Tf := conv{(0, 0), (|, 0), (0, |)} and Tg := 
conv{(l, 1), (|, 1), (1, |)} are visualized. The right-hand sides of the primal (1) and dual 
problem (3) are 


f{v) 



resp. 



dv 

dxi 


dx. 


This corresponds to fi = 0, = (xr/,0), gi = 0, = {XTg,0), where for w C 

denotes the characteristic function, i.e., Xuj{^) = 1 for x G (u and Xa)(x) = 0 for x G 
Figure 1 also shows some approximations of the primal and dual solution. The primal 
solution u has a line singularity along conv{(|, 0), (0, |)}, while the dual solution z has a 
line singularity along conv{(|, 1), (1, |)}. At the intersection of the lines with dQ, there 
are point singularities. 

Figure 2 (left) shows the typical convergence behavior for the estimators r]u,i and rjz/, 
the estimator product riu/riz,i, and the goal error \g{u) — g{U()\, where we used Algo¬ 
rithm A-C with 6 = 0.5. Similar results are obtained for other choices of 0 < d < 1 (not 
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Figure 2. Example from Section 4.5: Over the numbers of elements 
we plot the estimators rju/ and the estimator product rju/rjz/, as well 
as the goal error \g{u) — g{Ue)\ as output of Algorithm A-C with 0 = 0.5 
(left) resp. the estimator product for various 0 G {0.1,..., 0.9} as well as 
for 0 = 1.0 which corresponds to uniform refinement. 

displayed). The estimator product riu,iVz,i shows the optimal convergence rate of 0{N~^) 
as predicted by theory for p = 3 in 2D. 

Figure 2 (right) shows that all Algorithms A-C yield the optimal rate of convergence 
for a large range of values of 9 including 9 = 0.9. Uniform refinement corresponds 
to 9 = 1.0 and shows a suboptimal rate of 0{N~^) 

Figure 3 shows the numerical results for standard AFEM, which are based on adap¬ 
tivity for either the primal or the dual problem. In both cases, theory predicts optimal 
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Figure 3. Example from Section 4.5: Output of standard (non-goal- 
oriented) AFEM algorithms with 9 = 0.5, where adaptive mesh-refinement 
is steered only by the primal error estimator rju/ (top) resp. the dual error 
estimator rjz/ (bottom). Over the numbers of elements #7^, we plot the es¬ 
timators riu,i and the estimator product riu/rjz/, as well as the goal error 
\g{u)—g{U£)\ (left) resp. the estimator product for various d G {0.1,..., 0.9} 
as well as for 0 = 1.0 which corresponds to uniform refinement. 
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#738 = 1,022 #720 = 1,146 #720 = 1,094 #722 = 1,010 


AFEM (dual) 



Figure 4. Example from Section 4.5: Meshes generated by Algorithm A, 
B, and C as well as standard (non-goal-oriented) AFEM driven by the 
primal error estimator resp. the dual error estimator (from left to right) for 
9 = 0.5. 


convergence behavior for the related error estimator, at least if the adaptivity 

parameter 9 is sulficiently small; see, e.g., [18, 23, 15]. For all 9 G (0.1,..., 0.9}, we 
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Figure 5. Example from Section 4.5: For Algorithm A, B, and C as 
well as standard (non-goal-oriented) AFEM driven by the primal error 
estimator resp. the dual error estimator, we plot the cumulative num¬ 
ber of elements Acum := "^^75 necessary to reach a prescribed ac¬ 

curacy riu/riz/ < tol over 9 G {0.1,, 0.9} for p = 3 and tol = 10“^ (left) 
resp. p = 2 and tol = 10“^ (right). 


observe the optimal rate for the error estimator which drives the adaptive pro¬ 

cess. However, for the estimator product r]u,*Vz,* these strategies result in a suboptimal 
convergence rate 0 {N~‘^). 

In adaptive computations, the overall runtime depends on the entire history of adap¬ 
tively generated meshes. To better compare the various algorithms. Figure 5 shows the 
cumulative number of elements 

e 

(42) 

j=0 

which is necessary to reach a prescribed accuracy of riu,iVz,i < tol) versus 9 G {0.1, 0.9). 
The definition of Acum reflects the total amount of work in the complete adaptive process. 
Altogether, we compare five adaptive strategies: Besides Algorithm A-C, we consider 
standard AFEM based on the primal error estimator and standard AFEM based on the 
dual error estimator. For example, for a tolerance tol = 10“^ and p — 3, Figure 5 
(left) shows that Acum is smallest for Algorithm B-C for 9 = 0.8. Furthermore, we 
see that the goal-oriented algorithms A-C are superior to standard AFEM. Amongst 
the goal-oriented algorithms, because of having combined primal and dual refinement. 
Algorithm B-C are superior to Algorithm A, which only does one-sided refinement per 
iteration step. Furthermore, Algorithm B is at least competitive and sometimes even 
superior to Algorithm C. As visible in Figure 5 (right), for tol = 10“^ and p = 2, Acum is 
smallest for Algorithm B and 9 = 0.6. 

5. Goal-Oriented Adaptive FEM for Flux Evaluation 

5.1. Model problem. On a bounded Lipschitz domain 12 C with boundary 
r := dVt and for given A, /i G T^(12), and /2 G we aim to compute the weighted 

boundary flux 

(43a) g(u) := j(A\/u) -nAds, 
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where u is the solution to (32). For smooth u, g{u) can be rewritten as 


(43b) g{u) = / div{A'Vu)z dx + / AVu • S/z = a(u, z) — f{z) =: Nz{u) 

Jn Ja 

for all 2 ; G with z\r = A. Since the right-hand side is well-defined for u G 

this is a valid generalization of the fiux [28, Section 7]. Let 2: be the unique solution of 

the following inhomogeneous Dirichlet problem: 

z G H^{Q) with z|r = A such that a(v, z) = 0 for all v G Hq{Q). 

Then, it holds 

NAu) = -/( 2 ). 

5.2. Discretization. For a given regular triangulation % oi fl and a polynomial 

degree p > 1, let P^(7K) be defined as in Section 4.2. Consider := 

and Sq{%.) := fl Hl{Q). Let C* be the unique FEM solution of the homogeneous 

Dirichlet problem 

(44a) U,eSl{%) such that a(t/*, E*) =/(E*) for all E* G 5^(7;). 

Suppose that A G «S^(7o|r) := {E)|r : E) G S^{%)} belongs to the discrete trace space 
on the initial triangulation %. To approximate Nz{u) from (43), we let Z* be the unique 
FEM solution of 

(44b) Z* G S^{%) with Z*|r = A snch that o(Er, Z^) — 0 for all Et G 
and define 

(45) NzAU.) = -f{Z.). 

Lemma 25. There holds 

\Nz(u) — Nz,i,(Ui,)\ < CfluxlE — L*|Ei(q)||z — Z*|Ei(n), 
where 11^, denotes the FEM approximation of u from (2) and Cflux > 0 depends only on 

Proof. Since z — Z^, E there holds 

|iV,(u) - NzAU.)\ = \f{z) - f{Z,)\ = \f{z -Z,)\ = \a{u, z - Z,)\ 

= \a{u - U^,z- Z*)| < ||u- U^\\Hi(a)\\z - Z^\\hi{ci), 
where we used the definition of z and Z*. ■ 

5.3. Residual error estimator. The residnal error estimator for the primal problem 
remains the same as in (35), i.e., 

(46) := hlWCWU. - - div/2||i.,„ + ftrl|[(AVf/. + • nll|iz|sTnQ)- 

Since the inhomogeneons bonndary data satisfies A G «S^(7o|r) also the dnal estimator 
riz,* remains the same as in (36) with = 0 and g 2 = 0, i.e., 

(47) riz,i,{T)‘^ := -l- hrlli^AVZ* • 

Lemma 25 together with the reliability of for w G {u, z} (see, e.g., [3, Proposition 3] 
for the inhomogeneous Dirichlet problem for z) implies 

(48) \Nz{u)-Nz,,{U,)\<pu,.V,,*. 

Hence, the problem fits into the abstract framework of Section 2. We aim for optimal 
convergence of the right-hand side of (48). 
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Figure 6. Example from Section 5.5: The left figure shows the geome¬ 
try of the domain fl, the support of the primal Dirichlet data (blue), the 
direction of the primal convective field (blue arrow), the support of the 
dual Dirichlet data (red), and the direction of the dual convective field 
(red arrow). The right figure shows the initial triangulation To so that the 
inhomogeneous Dirichlet data belong to the discrete trace space <S^(7o|r)- 

5.4. Verification of axioms. With newest vertex bisection from [43] as mesh- 
refinement strategy, the assumptions of Section 3.2 are satisfied. It remains to verify 
the axioms (A1)-(A4), where dl^(7^ , %) := a{We - IF*, IF, - IF*)F2 ^ 

Theorem 26. Consider the model problem of Seetion 5.1. Then, the eonforming dis¬ 
cretization (44) of Section 5.2 with the residual error estimators (46)-(47) from See¬ 
tion 5.3 satisfies stability (Al), reduction (A2) withq^-g^ = 2“^/*^, quasi-orthogonality (A3), 
and discrete reliability (A4) with 7lu{%,%) — TZz{7i,%) = %\%. In particular, the Al¬ 
gorithms A-C are linearly eonvergent with optimal rates in the sense of Theorem 16, 17, 
20, and 21 for the upper bound in (48). 

Proof. For the primal problem, (A1)-(A4) follow as in Theorem 23. For the dual prob¬ 
lem, the axioms (Al)-(A2) follow from Theorem 23 since the estimator did not change. 
The discrete reliability (A4) is proved in [3] for general IF G iT^(r). In our particular 
situation, the proof simplifies vastly and shows even TZz{7i, %) = %\K. To see the quasi- 
orthogonality (A3), choose a discrete extension IF G S^{To) with IF|r = IF. Consider 
the solution G Sq{%) of 

a(F*, Z°) = -a(F*, IF) for all F* G S^{Z). 

Then, there holds Z* = -|- IF. Consequently, there holds dl^(7^^._^j ,Tej) — H-^q+i — 

||i,i(Q). Since is the solution to a homogeneous Dirichlet 
problem, the proof of (A3) follows analogously to that of Section 4. ■ 

5.5. Numerical experiment II: Flux-oriented adaptive FEM for convection- 

diffusion. We consider a numerical experiment similar to [36, Section 5.3] for some 
convection-diffusion problem in 2D. The goal of this experiment is to verify the optimal 
convergence of Algorithm A-C for the fiux quantity of interest (43) and, moreover, to 
illustrate this for a nonsymmetric second-order elliptic operator, which is covered by our 
theory. 

Let D = (0,1)^ C be the unit square. We set A = vl, with u > 0 the diffusion 
coefficient, b = {y,I — x), which is a rotating convective field around (|,0), and c = 0. 
According to div 6 = 0, it holds 

C = —uA -|- 6 • V and C^= —vA — 6 ■ V . 
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Algorithm A 



number of elements N = 






Figure 7. Example from Section 5.5: Over the numbers of elements 
we plot the estimators rju/ and rjz/, the estimator product as well as 

the goal error \Nz{u) — Nz/{U()\ as output of Algorithm A-C with 9 = 0.5 
(left) resp. the estimator product for various 9 G (0.1,... ,0.9} as well as 
for 9 = 1.0 which corresponds to uniform refinement. We consider p = 1 
and ly = 10“^ (right). 


We set f{v) =0 and consider non-homogeneous Dirichlet data on dfl for the primal 
problem, a pulse, defined by the continuous piecewise linear function 


UDiT{x,y) = < 6(| - x) 

[o 


l<x<l,y = 0 
l<x<l,y = 0 
otherwise. 
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10® 



- Algorithm B 

- Algorithm C 

10® I-^^^^^^^- 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

parameter 9 

Figure 8. Example from Section 5.5: For Algorithm A, B, and C as 
well as standard (non-goal-oriented) AFEM driven by the primal error 
estimator resp. the dual error estimator, we plot the cumulative num¬ 
ber of elements Acum := Xlj=o "^^75 necessary to reach a prescribed ac¬ 
curacy riu/riz/ < 10“^ over d G {0.1,..., 0.9} for p = 1 and u = 10“^. 




#7t = 337 



#Ti2 = 1,798 



#ri8 = 20,634 


Figure 9. Example from Section 5.5: Primal approximations Ui (top), 
dual approximations Zi (middle) and adaptively generated meshes 7} (bot¬ 
tom) for i G {6,12,18} (from left to right) as output of Algorithm B for 
9 = 0.6 and u = 10“^. Although we use a non-stabilized Galerkin scheme, 
initial oscillations in unresolved boundary layers are picked up immediately 
by the adaptive algorithm for both, the primal as the dual solution. 


Note that UDir trivially extends to some discrete function uoir G *5^(70) if % is cho¬ 
sen appropriately. Therefore, we can rewrite the problem into a homogeneous Dirichlet 
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#r22 = 28,839 



#r3i = 21,815 



#754 = 20,361 


Figure 10. Example from Section 5.5: To study the robustness of 
the goal-oriented algorithm with respect to the diffusion coefficient v G 
{10“^, 10“^, 10“^} (left, from top to bottom), we plot the estimators 
and the estimator product rju/rjz/, as well as the goal error 

\Nz{u) — N 2 /{Ui)\ as output of Algorithm B with 6 = 0.6 over the num¬ 
bers of elements (left). We show some related discrete meshes with 
> 20,000 elements (right). 

problem. To that end, write u = Uq + iiDir with Uq G Hq (fl) and solve 
a{uo,v) = f{v) — o(uDirW) ^11 V G Hl(yt). 

Note that the additional term on the right-hand side is of the form divA -|- A for some 

7o-element wise constant A and some A G L^(r2). A direct computation shows that 
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the weighted-residual error estimator with respect to uq coincides with Arguing as 
in the proof of Theorem 26, we see that the estimator satishes the axioms (Al)-(A4). 
Altogether, the problem thus fits in the frame of our analysis. 

The primal solution corresponds to the clockwise convection-diffusion of this pulse. 
We choose the boundary weight function A : dVi ^ E as a shifted version of the above 
pulse: 

(6{x-l) l<x<l,y = 0 

Hx,y) = <6(1 - x) l<x<l,y = 0 
I 0 otherwise. 

The dual solution corresponds to the counter-clockwise convection-diffusion of this pulse. 
For small u, the (primal and dual) pulses are transported from dQ into O and eventually 
back to dfl where a boundary layer develops. See Figure 6 (left) for an illustration of 
the supports of the primal and dual Dirichlet data, and the primal and dual convective 
fields. 

All discrete approximations are computed with lowest-order finite elements of de¬ 
gree p = 1. The uniform initial triangulation To is as shown in Figure 6 (right) ensures 
that the (primal and dual) Dirichlet data belong to the discrete trace space <S^(7o|r)- 

As shown in Figure 7, Algorithm A-C yield optimal convergence rates for the flux 
quantity of interest. For — 10“^ and a large range of values of 6* G {0.1,..., 0.9}, we 
observe the optimal convergence rate while uniform mesh-refinement appears to 

be slightly suboptimal. 

To compare the overall performance of the different algorithms, Figure 8 visualizes 
over different marking parameters 9 G {0.1,..., 0.9} the cumulative number of elements 
Acum which is necessary to reach a prescribed accuracy of riu,£Vz,i < 10“^ vs. the marking 
parameter d G {0.1,..., 0.9}; see (42) for the definition and interpretation of Acum- For 
Algorithms A-C, we observe that Acum is smallest for relatively large values 9 > 0.5, 
with Algorithm A being less efficient than Algorithm B and C. Overall, Algorithm B 
with 9 = 0.6 seems to be the best choice. 

Figure 9 shows several approximations and meshes obtained with Algorithm B. Be¬ 
cause u = 10“^ is relatively small, both the primal and the dual solution have significant 
boundary layers. These layers as well as the weak singularities coming from the kinks in 
the Dirichlet data are well captured by the adaptive algorithm. 

Figure 10 illustrates the effect of varying p G {10“^, 10“^, 10“®}. The optimal conver¬ 
gence rate of the estimator product is observed for the indicated values of p, however, 
the pre-asymptotic regime is longer for smaller values of p. This is to be expected, as the 
hidden constant in (48) depends on the reliability constant for the estimators, which in 
turn depends on p. 


6. Goal oriented BEM 


In this section, we extend ideas from [24] and prove that our abstract frame of convergence 
and optimality of goal-oriented adaptivity applies, in particular, to the BEM. 

6.1. Model problem. Let F C dfl denote some relatively open boundary part of 
the Lipschitz domain fl C E^, d = 2,3. Given F, A G iT^(r), we aim to compute the 
weighted boundary flux 




:= / Auds, 
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(49) 


where u solves the weakly-singular integral equation 


(50) 


Vu{x) 


G{x,y)u{y) dy = F{x) 


almost everywhere on F. 


Here, G : \ {0} —>• M denotes the Newton kernel 


G{x,y) 


1 

4:7T\x—y\ 


y\ for d = 2, 
for d = 3. 


The simple-layer operator extends to a linear and continuous operator V : H ^/^(F) —)■ 
where iF^/^(F) := {n|r : v G is the trace space of and iF“^/^(F) 

denotes its dual space with respect to the extended L^-scalar product; see, e.g., [35, 30, 41] 
for the mapping properties of V and the functional analytic setting. For d = 3 as well as 
supposed that diam(r2) < 1 for d = 2, the induced bilinear form 

a{u, v) := (Vu , v) := j {Vu){x)v{x) dx for u,v E X := 

is continuous, symmetric, and id“^/^(F)-elliptic. In particular, |||n||p := a{v,v) defines an 
equivalent norm on id“^/^(F). Moreover, the problem fits in the frame of Section 1.2. 
More precisely and according to the Hahn-Banach theorem, (50) is equivalent to (1), 
where the right-hand side of (1) reads f{v) := Jj^Fvdx. Moreover, the goal functional 
from (49) satisfies g E id“^/^(F)* = id^/^(F), where the integral is understood as the 
duality pairing between id“^/^(F) and its dual id^/^(F). 


6.2. Discretization. Let 71^ be a regular triangulation of F into affine line segments 
for d = 2 resp. flat surface triangles for d = 3. For each element T E 71, let ■jt '■ dlef ^ T 
be an affine bijection, where the reference element is Tref = [0,1] for d = 2 resp. T^ef = 
conv{(0, 0), (0,1), (1, 0)} for d = 3. For some polynomial degree p > 1, define 

T* := V^{%) := {K : F ^ E : K o 7 t G ^^(dlef) for all T E %}, 


where 7^^(Tref) := {q E L‘^{Tj.ei) : q is polynomial of degree < p on T^ef}. Let E 

T* be the unique BEM solutions of (2) resp. (5), i.e., 

(51a) G* G VP{%) such that a(t/*, E*) = /(E*) for all E* G V^iZ), 

(51b) Z,eVZ%) such that o(I4,Z*)-p(I/*) for all E* G 1P^(7;). 


6.3. Residual error estimator. The residual error estimators from [16] for the 
discrete primal problem (2) and the discrete dual problem (5) read 

(52) g^ZT?:=hT\\V{VZ-F)\\l,^^^ and := hrWViVZ, - A)\\l,^^y 

The error estimators satisfy reliability (6); see, e.g., [16]. The abstract analysis of Sec¬ 
tion 1.2 thus results in 

(53) \g{u) - g{Z) \ < 

and we aim for optimal convergence of the right-hand side. 

6.4. Verification of axioms. With 2D newest vertex bisection from [43] for d = 3 
resp. the extended ID bisection from [2] for d = 2 as mesh-refinement strategy, the 
assumptions of Section 3.2 are satisfied. It remains to verify the axioms (A1)-(A4), 
where dl^(7^,r*) := |||IF, - IF*||| ~ ||IF, - 
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parameter domain 0 < s < 2 x 


Figure 11. Example from Section 6.5 with conforming weight: Domain 0 
with initial triangulation To (upper left), exact solution u plotted over the 
boundary (upper right), exact primal and dual solution plotted over the 
arc-length, where s = 1 corresponds to the reentrant corner and s — 0.25 
corresponds to Zq (lower left), and adaptive mesh with #72o = 279 elements 
generated by Algorithm B with 6 = 0.5 (lower right). 


Theorem 27. Consider the model problem of Seetion 6.1. Then, the eonforming dis¬ 
cretization (51) of Section 6.2 with the residual error estimators (52) satisfies stabil¬ 
ity (Al), reduction (A2) with gred = 2“^/^'^“^^ quasi-orthogonality (AS), and discrete 
reliability (A4) with TZyj(Ti,lfi) = {T e 7r : T n IJC^\'^) 7^ consists 

of the refined elements plus one additional layer. In particular, the Algorithms A-C are 
linearly convergent with optimal rates in the sense of Theorem 16, 17, 20, and 21 for the 
upper bound in (53). 

Proof. The assumptions (A1)-(A2) and (A4) are proved in [24, Proposition 4.2, Propo¬ 
sition 5.3] for the lowest-order case. The general case is proved in [21]. The quasi¬ 
orthogonality (AS) follows from symmetry of a(-, •) and (A4); see Remark 11. ■ 

6.5. Numerical experiment with conforming weight function. Let P C 

denote the L-shaped domain shown in Figure 11 which satisfies diam(P) = l/\/2. On 
the boundary F := dTt, consider the function (f{x) := cos(2q;/3) for polar coordinates 

r{x),a{x) with origin (0,0). Consider the model problem (50) with 

F-.= (;c +1/2)0, 
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Algorithm A 




Algorithm B 




Algorithm C 




Figure 12. Example from Section 6.5 with conforming weight: Over the 
number of elements we plot the estimators and the goal error \g{u) — 
g{Ui)\ as output of Algorithms A-C for 9 = 0.5 (left) resp. the goal error 
\g{u) — g{U()\ for various 9 G (0.1,... ,0.9} as well as for 9 = 1.0 which 
corresponds to uniform rehnement (right). 


where /C : i/^/^+®(r) —)■ i7^/^+®(r), for all —1/2 < s < 1/2, denotes the double-layer 
potential, which is formally defined as {uy denotes the outer unit normal on F at y) 

IC(t>{x) := ^ />(?/) dy. 
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Figure 13. Example from Section 6.5 with conforming weight: Over dif¬ 
ferent values of 6*, we plot the cumulative number of elements iVcum := 
necessary for Algorithms A-C to achieve an error accuracy 
\giu)-g{U()\ < 10“®. 

For these choices, it is known [30, 35, 41] that (50) is equivalent to the Laplace-Dirichlet 
problem 

AP = 0 in fl subject to Dirichlet boundary conditions P — 0 on F, 
and the exact solution of (50) is the normal derivative of P, 

u(x) = dn^P{x) for all a: e F. 

We define the initial mesh % as shown in Figure 11. As weight function A G <S^(7o), 
we consider the hat function defined by A{zo) = 1 and A{z) = 0 for all other nodes 2 ; of 
To (the node Zq is indicated in Figure 11). 

For the lowest-order case p = 0 and 0 = 0.5 in Algorithm A-C, Figure 12 shows the 
convergence rates of the error estimators r 7 „, their product rjuflz) and the error in the 
goal functional \g{u) — g{U()\. Moreover, we compare the convergence rate of the error 
in the goal functional for different values oi 9 G {0.1,... ,0.9}. For either choice of 9 
and all adaptive algorithms, we observe the optimal convergence rate (#7})“^^^ for the 
respective error estimators as well as (#7})“^ for the error in the goal functional. 

For different values of 0 G (0.1,... ,0.9}, Figure 13 plots the cumulative number of 
elements Acum := necessary to reach a given error tolerance 10“®. We observe 

that for all three algorithms a large 9 ~ 0.8 seems to be optimal. Moreover, Algorithms B- 
C show comparable performance which is clearly superior to that of Algorithm A in the 
whole range of 9. 

6.6. Numerical experiment with non-conforming weight function. We consider 
the same setting as in Section 6.5, with the only difference that A is the characteristic 
function of Fa ^ F, i.e., A(x) = 1 on Fa and A(a:) = 0 on F \ Fa. We choose Fa ^ F as 
the part of F which is marked in red in Figure 11. This implies that the goal functional 
takes the form 



Note that A ^ P^/^(F) D P^(F), but only A G P^/^“®(F) for all £ > 0. In particular, 
g ^ i7“^/^(F)* = i7^/^(F). Consequently, this example is not covered by the theory of 
the previous sections. This is also reflected by the numerical results, if the adaptive 
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Algorithm A 


□ □ 



Figure 14. Example from Section 6.6 with non-conforming weight: Coun¬ 
terexample to show that rescaling is necessary. We plot the output of Al¬ 
gorithm A for 6* = 0.5 without rescaling of the estimators, i.e., rju = and 
Vz = vti with e = 0. We do not observe convergence at all. 



Figure 15. Example from Section 6.6 with non-conforming weight: Piece- 
wise constant approximation of primal and dual solution plotted over the 
arc-length (left), and adaptive mesh with #725 = 299 elements generated 
by Algorithm B for 0 = 0.5 and with the rescaled estimators (54) with 
£ = 0.3 (right). 


algorithms are naively employed; see Figure 14, where we do not observe convergence at 
all. 

To account for the fact that A ^ H^{T), we approximate A in each adaptive step by the 
continuous function Ai G S^{Te) defined by Ai{z) = 1 for all nodes z oi Te with 2 : G Ea 
and Ai{z) = 0 for all other nodes. Convergence A^ —?• A is assured by marking of the 
two elements where Ai is not constant in each adaptive step. Since there clearly holds 
IAr||_ffi/ 2 (r) —>■ 00 as £ —>■ 00 , we need to rescale the error estimators for the primal and 
the dual problem, respectively. Given e > 0, define 

(54) and y,(T)^=/.y'||V(VZ, - 

Since a thorough analysis is beyond the scope of this paper, we only provide a heuristic 
motiviation for this rescaling: With A^ ~ A, the error in the goal functional is estimated 
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Algorithm C 




Figure 16. Example from Section 6.6 with non-conforming weight: Over 
the number of elements we plot the estimators and the goal error 

\g{u) — g{U() \ as output of Algorithms A-C for 6 = 0.5 (left) resp. the goal 
error \g{u) — g{U()\ for various 6 G {0.1,... ,0.9} as well as for 9 = 1.0 
which corresponds to uniform refinement (right), where we use the rescaled 
estimators (54) with e = 0.3. 


by 


\g{u)-g{Ui)\ ^ \ jj,u 


Ui){Ai VZ£)ds\ < \\u ||if-i/ 2 +£(r) 11^^ 11 if 1 / 2 -^ (r) 


31 









Figure 17. Example from Section 6.6 with non-conforming weight: Over 
different values of 6*, we plot the cumulative number of elements A^cum := 
Ylk=Q necessary for Algorithms A-C to achieve an error accuracy 
\g{u) — g{Ui)\ < 10“®. We use the rescaled estimators (54) with e = 0.3. 



Figure 18. Example from Section 6.6 with non-conforming weight: We 
plot the cumulative number of elements Acum necessary to reach the toler¬ 
ance gu,£Vz,i < 10“^ for different values of 6* G {0.1,... ,0.9} and different 
scaling parameters e G {0.1,... ,0.4}. We observe that e ~ 0.3 gives the 
best performance. 


Since sup^gj^ ll^^llRi/ 2 -e(r) < oo, the last estimate is even rigorous and follows from ap¬ 
propriate Poincare inequalities; see, e.g., [14, 16]. 

For £ = 0.3, 9 = 0.5, and lowest order BEM p = 0, Eigure 16 shows the convergence 
rates of the error estimators 77 ^, their product rjuVz, and the error in the goal functional 
\g{u) — g{U£)\. Moreover, we compare the convergence rates of the error in the goal 
functional for different values of 0 G {0.1,..., 0.9}. Except for 0 = 0.9 and Algorithm C, 
we observe for either choice of 9 and all adaptive algorithms the optimal convergence rate 
(#” 7 })“^ for the error in the goal functional as well as the estimator product. 
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Figure 17 plots the cumulative number of elements A^cum = Ylk=o"f^'^k necessary to 
reach a given error tolerance 10“® for different values of 9 e {0.1,. ..,0.9}. We observe 
that for Algorithms A-C a large 9 ~ 0.8 seems to be optimal, whereas Algorithm B shows 
optimal behavior for 0.3 < 9 < 0.7. Overall, Algorithm B seems to be the best choice in 
this experiment. 


7. Conclusions & Open Questions 

7.1. Analytical results. We have derived an abstract framework to prove conver¬ 
gence with optimal algebraic rates for goal-oriented adaptivity for finite element methods 
and boundary element methods. While the analysis of prior works [7, 37] was tailored 
to the Poisson model problem resp. symmetric boundary integral formulations [20], our 
approach which is inspired by [15], is a priori independent of the model problems and 
covers general linear second-order elliptic PDFs and fixed order elements in the frame 
of the Lax-Milgram lemma. Following [18], our argument avoids the discrete efficiency 
and hence the interior node property of the mesh-refinement required in [7, 37]. Fol¬ 
lowing [23], our argument uses the concept of a general quasi-orthogonality which allows 
to work beyond symmetric problems and, unlike [36, 29], to avoid any assumption on 
the initial mesh %. As firstly observed in [3] and later used in [23, 15], the convergence 
and quasi-optimality analysis relies essentially only on reliability of the error estimator 
(see axioms (A1)-(A4)), while efficiency is only used to characterize the estimator-based 
approximation classes in terms of the so-called total error, i.e., error plus data oscillations 
(see Lemma 22). In addition to the algorithm from [37] (Algorithm A), we gave a thor¬ 
ough analysis for the algorithm from [7] (Algorithm C) without additional assumptions 
on the given data. Moreover, we proposed a variant of the algorithms from [37] and [29] 
(Algorithm B). All three algorithms are proved to be linearly convergent with optimal 
algebraic rates (see Theorem 16, 17, 20, 21), where theory guarantees linear convergence 
for all marking parameters 0 < 9 < 1, while optimal convergence rates are qualitatively 
guaranteed for 0 < d < d* (Algorithm A-B) resp. 0 < d < 0*/2 (Algorithm C) for some 
a priori bound 0 < d* < 1 which depends on the given problem. 

7.2. Empirical results. To underline our analysis, we considered three different 

problems: First (Section 4.5), we computed an example from [37] which considers finite 
elements for the Poisson model problem with some right-hand side / = div/ and goal 
function g = divgf for some piecewise constant vector fields Essentially 

for all choices of adaptivity parameters 0.1 < 9 < 0.9, we observed optimal convergence 
behavior of the goal-oriented adaptive algorithms, while standard adaptivity leads to a 
reduced order of convergence. Second (Section 5.5), we modified an example from [36] 
with a non-symmetric operator, where the goal is the evaluation of the flux for some finite 
element computation. All goal-oriented adaptive algorithms are robust with respect to 
the choice of 0.1 < d < 0.9. Finally (Section 6.5), we considered an example in the 
frame of the boundary element method, where the goal was some local flux evaluation. 
Again, all goal-oriented adaptive algorithms are robust with respect to the choice of 
0.1 < 9 < 0.9 and lead to optimal convergence behavior. Throughout, our observation 
was that the new algorithm (Algorithm B) leads to the best results with respect to the 
cumulative sum of elements (42) which seems to be an appropriate measure for the overall 
computational performance to reach a prescribed accuracy. Although we did not observe 
that Algorithm C leads to suboptimal convergence rates for large 9, where Algorithm A 
and B still are optimal, we note that this has been observed in [20] for the point evaluation 
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in boundary element computations which is a linear and continuous functional (and hence 
an advantage) of boundary integral formulations. 

7.3. Extensions & open questions. First, following the work of Mommer & Steven¬ 
son [37, Section 7.1], it is possible to use the extraction framework to apply our conver¬ 
gence and quasi-optimality results to compute point values. Second, arguing along the 
lines of [23], we think that it is possible to include (at least certain) nonlinear goal func¬ 
tional and nonlinear PDFs based on strongly monotone operators. As in [23], we note 
that the proof of stability (Al) and reduction (A2) might be challenging for higher-order 
elements p > 2, since even optimality results for standard AFEM for nonlinear problems 
are restricted to the lowest-order case p = 1; see, e.g., [10, 15, 17, 23, 27]. Finally and for 
the ease of presentation, we focussed on (homogeneous) Dirichlet conditions throughout 
our experiments. We note that the extension to mixed Dirichlet-Nenmann-Robin bound¬ 
ary conditions is easily possible; see [3, 15, 25] in the frame of standard AFEM. However, 
we note that our analysis currently requires that the Dirichlet data belongs to the coars¬ 
est trace space 5^(7o|r)- The main reason is that our analysis uses that the difference of 
solution and EEM approximation, i.e., u — Ui for the primal problem resp. z — for the 
dual problem, is an admissible test function. The latter fails for general inhomogeneous 
Dirichlet conditions. We believe that the rigorous analysis of this problem is beyond the 
current work and requires further ideas beyond those of standard AEEM [3, 15, 25]. 
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